Pade approximants for the ground-state energy of closed-shell quantum dots 
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Analytic approximations to the ground-state energy of 
closed-shell quantum dots (number of electrons from 2 to 
210) are presented in the form of two-point Pade approxi- 
mants. These Fade approximants are constructed from the 
small- and large-density limits of the energy. We estimated 
that the maximum error, reached for intermediate densities, 
is less than < 3 %. Within the present approximation the 
ground-state is found to be unpolarized. 

FACS numbers: 73.20.Dx, 73.20.-r, 71.45. Gm 



Quantum dots are artificial atoms which have been a 
subject of inteuse theoretical and experimental research 
in recent years.u They offer the very interesting possibil- 
ity of varying the parameters characterising the dot in 
a broad range, and also of fixing at will the number of 
electrons, N , confined in the dot. External magnetic and 
electric fields can be added to generate new static and 
dynamical effects. 

For very strong confinement (high electron densities) 
the one-particle Hamiltonian dominates over the coulomb 
repulsion, which can thus be treated as a perturbation. 
On the other hand, when the confinement is soft and the 
electron density becomes small, correlation effects start 
to play an important role, leading at very srriall densities 
to crystallisation,El as envisaged by Wigner.tl Signals of 
the crystal phase can be seen even in the j£vel structure 
of the N — 2 - A systems in one dimension,cl or in the "ge- 
ometry" of the wave functions (the spatial distribution of 
the electron probability) of the three-electron system in 
two dimensions Jj 

In the present paper, we exploit the fact that the 
ground-state energy can be computed analytically in 
both the large- and small-density limits, from which we 
construct an estimate for the energy for arbitrary elec- 
tron density by means of a two-point Fade approximant. 
The idea was applied Ions, ago to the three-dimensional 
(infinite) electron system.l3 We have used it recently to 
compute the low.- lying energy levels of < 5 electrons in 
a quantum dot,Ll and is extended here to dots containing 
as many as 210 electrons. 

One may ask why an interpolant should be a good ap- 
proximation to the ground-state energy. The reason is 
that the energy is a very smooth function of the density. 
For small systems, it was noticed in Ref. ^ that the low- 



density picture is valid up to surprisingly high densities, 
and in Ref. ^ it was stressed that the regions of conver- 
gence of both limiting expansions for the energy overlap. 
Below, we will argue that these regions also overlap for 
larger TV- values. 

We study quantum dots consisting of N electrons mov- 
ing in two dimensions under the action of a parabolic po- 
tential. The Hamiltonian describing the system may be 
written in dimcnsionless form as 



H 

hujo 



(1) 



where ujq is the dot frequency, (3^ = \J {-^-^) / (huJo) is a 
measure for the strength of the effective electron-electron 
interaction, fi is the electron effective mass, and k the di- 
electric constant of the material the electrons are moving 
in. We use the oscillator's natural units, i.e. y^h/ (/iWg) 
for length, hivo for energy, etc. Notice that /3 is the only 
parameter entering the hamiltonian h. By varying /3 we 
modify the "density" oLlie system. 

Following references,Q'u we will construct two-point 
Fade approximants from the series expansions of the en- 
ergy, e, in the /3 — > and /3 oo limits 
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in which bo is the energy of N free electrons in the 
quadratic potential, 63 4'olX]i<j ^ ^jl^^l^o >, 
ao is the classical (harmonic plus coulomb) energy of the 
iV-electron system, 02 = is the zero-point energy 

of the classical cluster, i.e. the uji are the small-oscillation 
frequencies around the equilibrium configuration. f3 —^ 
will be called the oscillator limit, and f3 ^ 00 the Wigner 
limit indicating that a finite- analog of the Wigner crys- 
tal (a Wigner "cluster") is formed in this limit. 

We computed the coefficients bo, 63, ao and 02 for dots 
with up to 210 electrons. Higher coefficients could, in 
principle, be obtained although the calculations become 
rather involved. From these coefficients we may construct 
the following approximants interpolating between the ex- 
pansions for /3 — » and /3 — > 00 



P3,2{l3) = bo + ao0Ul 
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where 92 = ao/(6o - ^2), 9i = aoq2/h, and 93 = (aogi - 
^3) / (&o — 02)- ^3.2 and P4_3 are the first two terms of the 
second-diagonal Pade sequence {Ps+i^s(/3)}. The labels 
s, t in Pg^t mean that we are taking s + 1 coefficients from 
the (3^0 series and t + 1 coefficients from the /3 — > 00 
series to construct the approximant. We will use (4) and 
(5) to estimate the ground-state energy of the iV-electron 
dot. 

In what follows, we restrict the analysis to closed-shell 
quantum dots, i.e. systems in which the first Nsheii os- 
cillator shells are filled for /? = which will simplify our 
calculation significantly. The number of electrons is thus 
restricted to the following values: = N shell {N shell + ^)- 
The energy in the /3 = limit is bo = N{2Nsheii + l)/3. 
For these systems, there is only one level starting from 
e = 60 at /3 = 0. Its angular momentum and spin are, 
respectively, L — and 5 = 0. This will be the ground- 
state for small (3. On the other hand, for very large values 
of (3 the leading contributions to the energy (i.e. the co- 
efficients ao and 02) are independent of the spin state of 
the system. Thus, we may undoubtedly construct the 
small-/? and large-/? expansions for the lowest rotational 
invariant (L = 0), unpolarised {S = 0) state. In the rest 
of the paper, this will be called the "ground" state, al- 
though it is possible that higher spirustates become lower 
in energy for intermediate densities.BllJ 

Let us compute the coefficient 63. The wave function 
for /3 = is given by the Slater determinant 



'^i(l)x+(l) 
0i(l)X-(l) 

0JV/2(1)X+(1) 
0JV/2(1)X-(1) 
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where we numbered sequentially the harmonic wave func- 
tions (orbitals) 0i, . . . , 4>n/2\ and x+^ X- are the spin-up 
and spin-down functions. The explicit form of the one- 
electron orbitals is 



(7) 



where Cfej/| = yJkl/l-K ik + |Z|)!], and for the composed 
index n — {k, I) we defined the following sequential order: 
{(0,0)}, {(0,1), (0,-1)}, {(1,0), (0,2), (0,-2)}, . . . , etc, i.e. 
increasing |^| inside a shell with positive-Z orbitals are 
taken first. The energy corresponding to <j)k,i is Afe^; = 
1 -f 2fc + |Z|. In terms of these energies we have 60 — 
X^n '^k.i- The coefficient 63 can be written as 

63 = ^/(ni,ni) -1-4 ^ /(ni,n2)-2 ^ 7(711,^2), 



ni<n2 



ni<n2 



where 
/(n^,n2) = 

J(ni,rt2) = 
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l'^«i(ri)|^|(?!>«2(^2) 
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C iri)K, {r2)<Pn. {ri)^n, (r2), (10) 
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Both I{ni,n2) and J{ni,n2) can be reduced to sums of 
products of Gamma functions. 

The numerical results for the coefficients 63 are pre- 
sented in Table |. Notice that a ainiple Thomas-Fermi 
estimate for the present problemllil shows the scaling 
property: eTF{N,(3) = N^/^eTpi^, N^'^^ P) and, con- 
sequently, the 63 coefficient should exhibit a A^^/* be- 
haviour in the large- TV limit which is satisfied by the re- 
sults shown in Table |. 

The classical energies of the iV-electron clusters, i.e. 
the coefficients ao , were computed by vaesM^ pf a Monte 
Carlo technique with Newton optimisationJi3'E2l The val- 
ues are also shown in Table |. For large N, the 
Thomas-Fermi approximation predicts the dependence 
Co ~ 1.062 N^/^ . We may obtain another estimate for oq 
from the classical energy of an hexagonal lattice struc- 
ture, yielding ao « 0.968 N'"/^. In Table |, it is shown 
that the coefficients ao of the finite systems follow the 
law when N > 30. 

Once the configuration of minimum energy of the clas- 
sical problem is found, the small-oscillatioia-|problem is 
solved in order to compute the frequenciesEj and, con- 
sequently, to find a2. The results are shown in the last 
column of Table | A very rough estimate of the zero- 
point energy is obtained from the dispersion relation of 
the longitudinal mode of an hexagonal lattice in the (1,0) 
direction, yielding a2 w 0.78 /V^/"*. The resufis in Table | 
make evident the N^^^ law for N > 30. 

An evident conclusion that can be extracted from Ta- 
ble | is that the properties of quantum dots with N > 100 
are dictated by the — > 00 asymptotics. With other 
words, the ground-state energy for such a dot should be 
accurately described by a P^^^ approximant in which 60 ~ 

|iV3/2, 63 « 0.7 ^ i Q ^5/3 g^^^ ^2 ^ Q g 

In this sense, these dots are more "infinite" than "meso- 
scopic" , and their properties will be very well reproduced 
by statistical theories. 

Another very general property following from the re- 
sults shown in Table | is the fact that the regions of 
convergence of the expansions (2) and (3) overlap. As 
an example, we show in Fig. y the P4_3 approximant 
(solid curve) along with the small-/? and large-/? ex- 
pansions for the energy (dashed curves) for the quan- 
tum dot with 42 electrons. The asymptotic curves al- 
most intersect one another at /? « 0.8. This smooth- 
ness of the dependence of e on /? guarantees that an in- 
terpolant will be a very good approximation to the ac- 
tual energy at any intermediate value of /?. The same 
property was shown to holdQ for the ground-state and 
the excited levels of small dots with N < 5. For large 
dots, this seems to be a general property which follows 
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from the scaling laws obeyed by the coefRcients. Indeed, 
we can obtain a naive estimate of the radii of conver- 
gence of the series (2) and (3) by comparing consecutive 
terms. The small-/? expansion is expected to work when 
/3 < /3o (60/63)^/^ « A^"i/i2^ ^jjjjg ^jjg large-/? expan- 
sion will be valid for (3 > (3^ ^ (02/00)^/^ « 7V-5/24^ 
Because /3o > Poo, the regions of convergence will over- 
lap. 

We show in Fig. |^ the "convergence" of the Fade se- 
quence, i.e. the relative difference between the P4_3 and 
the P3^2 approximants for the 42-electron system. This 
figure suggests that the relative error of the P^^^ approx- 
imant should not exceed 3 % over the entire interval 
< /? < 00. The ijxattern is similar to the one encoun- 
tered in small dotsj3 and therefore we expect this to be 
an estimate of the actual error. 

Finally, we turn ourselves to the question of the ab- 
solute ground-state of the system. In the present ap- 
proximation, the energy for large f3 is independent of the 
polarisation state of the system. The leading approxima- 
tion to the energy is again given by Eq. (2), where Oq 
and 02 are shown in Table ^ For intermediate (3 values, 
however, spatially antisymmetric states may be favoured 
as they minimise the Coulomb repulsion. We may easily 
compute the coefRcients &o and &3 at least for the lowest 
antisymmetric state at /3 << 1. This state is built up by 
placing one electron per orbital. If there are free orbitals 
in the last shell, the occupancy of the orbitals will lead 
to a maximal angular momentum (Hund's rule) . For the 
coefficient 63 we have now the expression 



To increase the usefulness of the present results we 
were able to fit the coefficients in Table I as follows 



&3 = 



ni<n2 



{/(ni,n2) - J(ni,n2)}, 



(11) 



where the sum runs over the occupied states. 

In Table ||, the angular momentum of such antisym- 
metric states and the coefficients 60 and 63 are given as 
function of N . Fade approximants are constructed in the 
same way as for unpolarised states. 

We show in Fig. ^ a comparison between the estimated 
energies of the polarised and unpolarised states for the 
dot with 42 electrons. Both curves become very close 
for /3 > 1, but the unpolarised state is always the low- 
est state for any /3. We found that this is the case for 
all N-values shown in Table ^ which implies that in 
the present approximation the unpolarised state corre- 
sponds to the absolute ground-state for closed-shell quan- 
tum dots. This agrees with recent quantum mechanical 
density functional calculations .llj 

In conclusion, we have obtained Fade approximants to 
the ground-state energy of N electrons (2 < < 210) 
which move in two dimensions and are confined by a 
parabolic potential. The approximants are asymptoti- 
cally exact in both the high- and low-density limits. The 
maximum relative error for intermediate densities is esti- 
mated to be only a few percent of the actual energy. This 
is also a typical accuracy reached in density.- fimctional 
calculations when the correlations are strong.olij 



ao/N^'^ = 1.062 - 0.875/A^^/2 - 0.185/A^, 
a2/iV5/4 ^ 0.573 + QAlb/N^'"^ - 0.160/A^, 
6o/iV3/2 ^ 2/3 + 0.083/A^, 
63/iV^/'^ = 0.731 - QAll/N^/'^ - 0.046/A^. 



(12) 

(13) 
(14) 
(15) 



For > 6 the fit to all the coefficients has a relative 
error which is smaller than 0.2%. 

In principle, the Pade-approximant method could 
equally well be applied to other dot configurations (i.e. 
non parabolic confinement, quasi- two-dimensional sys- 
tems) or t0|-Qt|her problems sucit,a&- for example, metal- 
lic clusters Jl3ii3 ions in traps ,t3~Ell etc. In these cases, 
more numerical work will be required to compute the 
coefficients &o, ^3: 10 and 02. 
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TABLE IL The angular momentum and the coefBcients 6o 
and ba for the lowest spatiaUy antisymmetric (spin polarised) 



state. 


N 


J 






2 


1 


1.06066 


0.186307 


6 





0.952579 


0.378007 


12 


6 


0.962250 


0.428734 


20 


5 


0.950329 


0.486891 


30 


12 


0.949386 


0.508812 


42 


18 


0.947863 


0.525187 


56 


10 


0.944959 


0.538148 


72 


36 


0.946083 


0.544753 


90 


12 


0.943998 


0.555149 



FIG. 1. Small- and large-/3 expansions (dashed curves) and 
the P4,3 approximant (solid curve) for the ground-state of the 
42-electron system. 



FIG. 2. Relative difference between the and the P^^z 
approximants for the 42-electron system. 



FIG. 3. Energy estimates for the polarised (dashed curve) 
and unpolarised (solid curve) states for a quantum dot with 
42 electrons. 



TABLE I. The coefficients bo, 63, ao and 02 for closed-shell 
quantum dots. 



N 


bo/iV3/2 




ao/iV^/^ 




2 


0.707107 


0.372611 


0.375 


0.784567 


6 


0.680414 


0.531250 


0.674147 


0.738995 


12 


0.673575 


0.591525 


0.793317 


0.698756 


20 


0.670820 


0.623655 


0.856344 


0.671474 


30 


0.669439 


0.643709 


0.895843 


0.653428 


42 


0.668648 


0.657442 


0.922532 


0.639116 


56 


0.668153 


0.667452 


0.941733 


0.633269 


72 


0.667823 


0.675071 


0.956422 


0.626248 


90 


0.667592 


0.681068 


0.967894 


0.620456 


110 


0.667424 


0.685913 


0.977103 


0.616908 


132 


0.667298 


0.689909 


0.984691 


0.613269 


156 


0.667201 


0.693262 


0.991046 


0.610032 


182 


0.667124 


0.696118 


0.996267 


0.607227 


210 


0.667063 


0.698572 


1.00107 


0.605239 
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Figure 2 
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